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ABSTRACT 

We study the dynamics of the chiral phase transition expected during the 
expansion of the quark-gluon plasma produced in a high energy hadron or heavy 
ion collision, using the 0(4) linear sigma model in the mean field approximation. 
Starting from an approximate equilibrium configuration at an initial proper time 
T in the disordered phase, we study the transition to the ordered broken symme- 
try phase as the system expands and cools. 

We give results for the proper time evolution of the effective pion mass and 
for the pion two point correlation function. We investigate the possibility of dis- 
oriented chiral condensate being formed during the expansion. In order to create 
large domains of disoriented chiral condensates low-momentum instabilities have 
to last for long enough periods of time. Our simulations show no instabilities for 
an initial thermal configuration. For the far-of-equilibrium cases studied, the in- 
stabilities are formed during the initial stages of the expansion survive for short 
proper times. For slow expansion rates even such configurations do not develop 
instabilities. 



*This work was done in collaboration with Fred Cooper, Alex Kovner, Emil Mottola and Juan Pablo 
Paz. This research was performed in part using the resources located at the Advanced Computing 
Laboratory of Los Alamos National Laboratory, Los Alamos, NM 87545. 



1. Model and Approximations 

The possibility of producing large correlated regions of quark condensate 
< QiQj > pointing along the wrong direction in isospin space was proposed 
to explain rare events where there is a deficit or excess of neutral pions observed 
in cosmic ray experiments'^. The idea that such disoriented chiral condensates 
(DCC's) can be formed has been the source of several experimental proposals^ in 
high energy collisions. It has been recognized that nonequilibrium dynamics can 
yield regions of DCC^~-^°. These authors introduced instabilities of low momenta 
modes by quenching the system, i.e., evolving a "hot" state under the dynamics 
of a zero temperature hamiltonian. 

There are clearly two important questions that the theoretical models should 
aim to answer. The first one is to determine whether during the evolution that 
follows the collision there are instabilities affecting the fluctuations, in which case 
there is a chance of the correlations growing. 

The second question is, assuming the instability exists, can the correlated 
domains grow large enough so that many pions can be emitted from each do- 
main making the detection of DCC's possible. Experimentally, one can look at 
the probability distribution of the fraction of neutral pions P(/) where / = 
n^ojin^o + n^±)- Assuming that all isospin orientations of tt are equally dis- 
tributed so that P(fi) = const, and using the equation P(f2)c?f2 = P(^f)df, one 
finds that for an ideal case where only one large domain is formed this distribution 
is given by P{f) = 1/2-^/, provided that a |7r*p, where / = cos^ 6. If the final 
particle distribution results from uncorrelated events, one would expect the neu- 
tral fraction to be distributed binomially around / = 1/3. Thus, for an ensemble 
with domains of different sizes the resulting distribution P(/) dwells somewhere 
in between binomial and domain-like distribution. To address these questions we 
employ the 0(4) linear sigma model, in the mean field approximations, which 
seems to have the essential attributes of being simple but realistic enough: it 
has the correct chiral symmetry properties and also appropriately describes the 
low energy phenomenology of pions. Moreover, the mean field approximation im- 
proves the classical approximation, because it takes into account quantum fluc- 
tuations that can not be neglected when the system cools down. In this model 
the mesons are organized in an 0(4) vector (^i, ^3) ^4) = {(^,'^i,'^2,'^3)- The 
effective action is derived by using the the large N expansion of a 0(iV) model, 
and is given by 

r[$, x] = / d'x{-^ J2 ^»(° + X)^. + ^ + \xv^ + Ha+ '-Nil InGo^}, (1) 
1=1 

where G^\x,y) = i[n + x{x)] 6\x - y), and x = A(E^=i " v') with iV = 4. 
Varying the effective action with respect to the mean fields, we derive their 
equations of motion 



(2) 



N 

X{x) = \{-v' + ^[^, \x) + Go{x, x)]). (3) 

i=l 

The functions Go{x,x) that appear in (3) are the coincidence limit of the prop- 
agators Go{x,y) that invert the operator Gq^ . 

By calculating two and four point function the parameters A, H and v were 
chosen to give the best fit to the physical observables /„-, and the phase 
shiftsi°. 

In Bjorken's Baked Alaska picture boost invariant initial conditions are fixed 
at an initial proper time tq (these initial configurations eventually lead to a flat 
rapidity distribution of pions as expected in the future RHIC experiments). Dy- 
namics then incorporates a cooling mechanism that may lead to low momentum 
instabilities (in contrast to a crude quench approximation that plants the insta- 
bilities from the beginning). In the boost invariant picture it is natural to use 
the following set of coordinates 

r^(^2-^2)l/^ ^^llog(i±^), x^, (4) 

z t — z 

where r is the proper time and rj is the spatial rapidity. The equations of motion 
reduce to 

T-^d^ rdr Ur) + X{r) Ur) = HSn, (5) 
N 

x{t) = \{-v' + \t)+ < ^\r,, x^,t) >]), (6) 

i=l 

(r-^drTdr -T-^dl-dl+x{T))v{ri,XuT) = Q, (7) 

where we have introduced an auxiliary field (p(^x,T), related to the Wightman 
function by: 

Go{ri^,x_i,T^,riy,y_i,Ty) =< ip{rix,x±,T^) ip{riy,y_i,Ty) > . (8) 

At equal proper times 

Go{ri^,x_i,riy,y_i;T) = Gc{rix,x±,riy,y_i;T), (9) 

where Gc{x,y) = jj^ Ylf=i < 4>i{^)4>i{y) > -4>i{^)4>i{y)- We expand the field ^p in 
an orthonormal basis 

(/?(?/, a;_L, r) = j dk{exp{ik ■ x)/fc(T) + h.c), (10) 

where A; • x = k^jT] + k±x±, dk = dkjjd^k±/[2'K)'^, and the mode functions /^(t) 
evolve according to: 

/, = -(^ + kl+ x(r) + ^ -ujlh. (11) 



A dot here denotes the derivative with respect to the proper time r. The Fourier 
transform of the expectation value < ((/3)^(a;,T) > can be expressed in terms of 
the mode functions fk and of the distribution functions 

Uk =< alak >, gk =< ako-k >, gl =< 44 > (^2) 

which entirely characterize the initial state of the quantum field. 

In order to solve equations (5), (6) and (11) as an initial value problem, 
we need to fix at tq the mean values (j)ik, 4>ik: ^fc) 9k and g^. This amounts to 
specifying five initial conditions at each momentum A;. * In the boost invariant 
case 4'ik{T) = 4'iki'^) = k ^ 0. 

The complexity of the processes we study suggests that our initial configu- 
ration is not a pure state but a mixed ensemble that should be described by a 
density matrix. In the mean field approximation the most general density matrix 
for each mode A; is a gaussian with five parameters associated with the above 
initial conditions and it can be written as 

P(.(t^ik,(t)ik) = ^2xGfc ^^P"i^'^'fe(^'fe ~ ^^'^^ ~ - 4:G^k^^k)\(t^^k " ^^^fc P 

- ^jfcp - l^^^i. - ^jfcp] - —-\(l)ik + cl>[k - 2^ifcP}, (13) 

0<^ik 

where 

Gjfc = < cl)ik cl>*k > -hk ^*k (14) 

Kik = < (t^zk (l>*k > -4>ik 4>zk (15) 

A'^ikGik = < {(l>%k,(l>*k] > -{4>%k,(l>ik]- (16) 

Pure density matrix satisfies the condition p = K — 4GS^ + 1/4G = (no 
4>i{k)(j)'^[k) mixed term). This subset corresponds to squeezed-coherent states-^-^. 
For S = and Gik = l/\/2wfc, the density matrix operator becomes p'^ = \c >< c\ 
where |c > is a coherent state. One can make a direct connection between our 
calculation and earlier work^~^, which used classical approximation to the same 
model. The equilibrium density matrix of a free theory is Gaussian and in the 
high T limit {Pu)k << 1) has the following representation: 

P=]lJ dhkd^^k exp{-|(|^,fc|' + {k^ + rr?)\^,k?)}p1j^{^,k, ^.J. (17) 

ik 

Therefore, solving the mean field equations with an initial thermal homogeneous 
density matrix is equivalent to solving the semi classical equations for a thermal 



*The formulation of Eqs. (5), (6) and (11) may give an impression, that we have to specify extra four 
initial conditions for /^(to) and fk{To)- This impression is false, since as it can be easily checked that 
the change of /^(to) and /^(to) if accompanied by an appropriate change of nj, and according to a 
Bogolyubov transformation leaves G^, and G^ invariant, and these quantities (together with 4'ik, 
<f)^^) unambigously specify the initial data. Therefore, we will fix the initial data fk{To) and /^(to) so 
that the vacuum state coincides with the ordinary Minkowsky vacuum, at least for high momentum. 



ensemble of inhomogeneous "classical" initial conditions. We want to stress the 
point that in the framework of the the mean field approximation one has to 
fix five initial conditions to describe the whole ensemble, instead of separately 
calculating the equations of motion for each member of the ensemble, as has been 
done in^~^. 

In the interacting case, p*''^ = e~^^ is not a gaussian any more, and therefore 
it is not contained in the set of mean field density matrices. However, one can 
approximate the thermal density matrix by replacing the quartic terms by using 
the thermal mass. We note that for a gaussian density matrix the requirement 
of isospin invariance leads immediately to a factorization form p = Yli Pi- Due 
to the factorizability of the density matrix, if one species i has a large average 

< iV* >, the whole density matrix will also have a large < N > and vice versa. 
This is also true for the variance < (iV — iV)^ >. The mean field approximation 
can only accommodate ensembles with either large fluctuations in both the to- 
tal number of pions and number of neutral pions or small fluctuations in both 
observables. Ensembles that have sharply defined total number of particles but 
large fluctuations in the number of neutral pions can not be described in this 
framework. 

2. Results 

In Fig. 1 we display the time evolution of the effective mass X: the mean den- 
sity of particles < n > = < dN/ drjdxj_ > and the relative fluctuation oi dN/ drjdxj_ 
starting at initial proper time tq = 1/m/c. We also plot the time evolution of 
the fourier transform of the correlator Gc{Arj,Axj_;T) at the planes |kj^| = 
and kr, = 0. The left row corresponds to an initial thermal configuration. In 
this case x does not become negative during the evolution so that no instabil- 
ities develop. The density < dN/ drjdxj_ > grows very mildly and so does its 
relative fluctuation. The fluctuation at large proper times is almost equivalent 
to its initial value. Therefore, observing a DCC for an initial thermal ensemble 
is unlikely. From the inverse width at half height of the profile of the Fourier 
transform we learn that the size of a domain in the beam direction is about one 
unit of space-time rapidity and about 1/m or less in the transverse direction. In 
the right row we generated instabilities by choosing an initial ensemble that is 
far from equilibrium for which we modified only the initial condition for a from 
to —1. We see that x becomes negative at early stages of the evolution and at 
this time an exponential growth of the low momenta modes occur. The density 

< dN/ drjdxj_ > triples itself by sucking energy from the a field and reaches a 
plateau after few tens of fermi/c. The relative fluctuation of this density also 
grows substantially and therefore the observation of DCC is more likely for this 
"peculiar" ensemble. This is also reflected in the structures appearing in the 
fourier transforms of Gc- More specifically we see that at the transverse direction 
there is an amplification of the amplitude at low momenta, which corresponds 
to a domain of about 2/m. However, in this case it is hard to define the size of 
the domain by naively measuring the inverse width due to multiple scales that 



emerge during the dynamics. For similar out of equilibrium ensembles we found 
no instabilities if we start the simulations at proper times greater than 2/m/c. 
This implies that if the expansion rate is too slow or if we start with lower energy 
densities even for out of equilibrium ensembles, DCC will be hardly produced. 
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Fig. 1. Proper time evolution of the x field, < dN/drjdxj_ >, the relative fiuctuation of 
dN/dr]dx± times the square root of the extent of the system (Accj^A?/)-^/^ and of the fourier 
transforms of the Green's function Gc at the planes oi p = |kj^| = and k^, = 0. The left 
row is for an initial thermal distribution where tc^ = l/(exp/3a;fc — 1), ct(to = 1) = ^eq{T = 
200MeV), 7f*(l) = = ^■(l) = 0. The right row has the same initial conditions as in the 

left one, but for o'(l) = —1. 



